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ABSTRACT 

We show that Type la supernovae (SNe la) are formed within both very young and old steUar 
populations, with observed rates that depend on the stellar mass and mean star-formation rates 
(SFRs) of their host galaxies. Models where the SN la rate depends solely on host galaxy stellar 
mass are ruled out with >99% confidence. Our analysis is based on 100 spectroscopically-confirmed 
SNe la, plus 24 photometrically-classified events, all from the Supernova Legacy Survey (SNLS) and 
distributed over 0.2<z<0.75. Using multi-band photometry, we estimate stellar masses and SFRs 
for the SN la host galaxies by fitting their broad-band spectral energy distributions with the galaxy 
spectral synthesis code, PEGASE.2. We show that the SN la rate per unit mass is proportional to the 
specific SFR of the parent galaxies ~ more vigorously star-forming galaxies host more SNe la per unit 
stellar mass, broadly equivalent to the trend of increasing SN la rate in later-type galaxies seen in the 
local universe. Following earlier suggestions for a simple "two-component" model approximating the 
SN la rate, we find bivariate linear dependencies of the SN la rate on both the stellar masses and the 
mean SFRs of the host systems. We find that the SN la rate can be well represented as the sum of 
5.3 ± 1.1 X 10~^^ SNe per year per unit stellar mass, and 3.9 ± 0.7 x 10"'' SNe per year per M© yr~^ 
of star formation. 

We also demonstrate a dependence of distant SN la light-curve shapes on star-formation in the host 
galaxy, similar to trends observed locally. Passive galaxies, with no star-formation, preferentially host 
faster-declining/dimmer SNe la, while slower-declining/brighter events are only found in systems with 
ongoing star-formation. We model the light-curve width distribution in star-forming galaxies as the 
sum of a young component, and an old component taken from the distribution in non-star- forming 
galaxies. Empirically understanding these relationships between SNe la and their environments will 
lead to future improvements in their use as cosmological candles. 
Subject headings: surveys - supernovae: general - distance scale - galaxies: evolution 



1. INTRODUCTION 

Type la supernovae (SNe la) represent cosmolo- 
gists' most direct probe of the cosmic expansion his- 
tory, yet an understanding of the composition of their 
progenitor systems has not yet been achieved (e.g. 
iHillebrandt fc Niemeveill200(t) . In principle, this uncer- 
tainty can be reduced and constraints placed on the na- 
ture of the progenitor if the typical explosion time-scale 
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of SNe la can be determined. This "delay-time" param- 
eterizes the distribution of times between the binary sys- 
tem formation and subsequent SN explosion following 
accretion of material from a secondar y co mpanion (e.g. 
iMadau. Delia Valle. fc PanagiallT998al see lGreggioll2005l 
for a review of current rate model formalizations). 

Various pieces of observational evidence have been used 
to place different constraints on the value of this delay 
time. The principal approach compares observed SN la 
rates with a predicted rate generated by convolving a de- 
lay functi on with an assumed cosmic star- formation his - 
tory fe.g. IMadau 'eranil998at iGal-Yam fc Maozll200l . 
Studies using this approach have determined a wide 
range of delay-t imes: 2±2-4Gyr ( Strolger et al.l 12004 
Imm. > 2 Gvr (IGal-Yam fc Ma.oH 120041) . and < 1 Gyr 
ijBarris fc Tonrvll2006(l . 

A second technique uses comparisons of the SN la rate 
in z < 1 galaxy clusters with the observed cluster iron 
content. If SNe la are assumed to represent the dom- 
inant source of iron in clusters, the low cluster SN la 
rate at low redshift implies much of the iron must have 
been produced from events at higher redshift - and hence 
suggests delay times of <2 Gyr (Maoz fc Gal- Yam 2004J. 
The third piece of evidence follows from the observa- 
tion that SNe la are substantially more common in star- 
form ing later-type galaxies t han in early-type systems 
Ce.g. lOemler fc TinslevllT979ll . The rate" per unit mass 
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is significantly higher both in later-type galaxies than in 
E/Sn svstems llva,Ti den Berghlll 99(11 iD^itTValle T.ivinI 
119941: iMannucci et a l.M2005l). and in bluer galaxies than 
in red galaxies ijMan nucci et al. 2005), with an enhance- 
ment of SNc la in ear ly-type galaxies that are radio- 
loud (Deli a Valle et al.ll2005.) . Furthermore, within star- 
forming galaxie s, SNe la are rarer in galaxy bulg es than 
in their discs ()Wang. Hoflich. fc Wheelal Il997|) . This 
suggests some dependence of the SN la rate on recent 
star- formation (and hence very short delay-times). Fi- 
nally, there is the observation that SNe la are common 
in old ev olved systems with litt l e recent star-formatio n 
activity llCaDDellaro et alJ 119991: IMannucci ei~an 120051) . 
suggesting some progenitors have delay-times of at least 
several Gyr; other studies of the star-formation histories 
of local SN la host galaxies claim a delay-time lower limit 
of -2 Gyr (Gallag her et al.ll200l . 

These contradictions can be resolved by removing 
the constraint of a single delay-time parameterizing 
all SN la explosions, and instead using a "two com- 
ponent" distribution (or even a more genera l func- 
tion) simila r to th e models proposed by IM annuc ci et al.1 
(|2005l 120061) and iScannanieco fc Bildstea (;200,5i. here- 
after SB05). Such models comprise a "prompt" (or 
small delay-time) SN la component, essentially depen- 
dent on recent star-formation, and an "old" (larger delay- 
time) component, dependent on the number of lower- 
mass stars. The total rate of SNe la is then a combi- 
nation of these different functions. These scenarios are 
able to resolve many of the observational contradictions 
described above (see discussion in SB05). 

These two-component models have implicit but impor- 
tant implications for the progenitors of SNe la which 
might impact their use as calibrated standard candles 
to derive cosmological parameters ('Ri ess et al.l Il998t 
[Perlmuttcr ct al. 1999; Asticr et al. 200$^ The possi- 
bility of subtle differences between SNe la from the two 
components, or a change in the relative fraction of the 
two types with redshift, are potentially concerning. As 
such it is vital to test and parameterize the SN la mod- 
els in as many ways as possible. One way to do this is 
to look at the rates and properties of SNe la in relation 
to the environments or galaxies in which they explode. 
Galaxies across the Hubble sequence provide an ideal lab- 
oratory for studying SNe la due the range of potential 
progenitor stellar populations that can be probed, from 
star-burst galaxies with dominant young populations of 
stars, through normal galaxies such as the Milky Way 
with a substantial fraction of evolved stellar mass, to the 
old, evolved elliptical galaxies in massive galaxy clusters 
which are essentially comprised of homogeneous old stel- 
lar populations. 

In this paper, we examine the properties of high- 
redshift SN la host galaxies, and calculate the frequency 
of occurrence of SNe la in galaxies of different type, ex- 
ploring the parameters governing the SN la rate. Our 
motivation is to measure the SN la rate as a function 
of the stellar mass and star-formation rates of the host 
galaxies in order to test the various predictions of the 
models described above. Performing these tests requires 
not only a large, homogeneous dataset of SNe la with 
well-understood detection efficiency characteristics, but 
additionally a multi-wavelength dataset for their host 
galaxy systems which can be used to constrain their 



stellar populations. The intermediate-redshift Supernova 
Legacy Survey (SNLS) satisfies both requirements. 

The SNLS uses repeat u*g'r'i'z' imaging of four 
square-degree fields to conduct a "rolling" high-redshift 
SN search. The repeat imaging allows not only the con- 
structio n of high-quality multi-band g'r'i' z' SN light- 
curves ()Astier et al.ll200(^ . but also extremely deep im- 
ages of the survey fields, and hence precise information 
on the broad-band spectral energy distributions (SEDs) 
of the SN host galaxies can be collected. Perhaps surpris- 
ingly, precision work on SNe la rates is most efficiently 
performed at these intermediate redshifts [z ~ 0.5); low- 
redshift SN la rates ar e still far less pr ecise than those 
at higher redshift (e.g. Neill et al. 2006j) due to the ho- 
mogeneity of high-redshift rolling SN searches. 

A plan of the paper follows. § El § El and § ^ detail 
the framework upon which our subsequent analyses in 
§ O and § are built. In § (21 we introduce the SNLS 
and discuss the construction of deep optical stacks of the 
survey field as well as the host galaxy identification and 
flux measurement. §0 describes the galaxy SED fitting, 
and §21 details the various incompleteness corrections we 
apply to our galaxy and SN la samples. Our analysis is 
contained in § [SI and §0 In § O we examine the rate of 
SNe la as a function of various host galaxy parameters, 
and in § examine the impact of environment on SN la 
light-curve shape parameters. We summarize in §[3 

Throughout the paper, we assume a cosmology of 
VIm = 0.3, f^A = 0.7, Ho = 70 kms"^ Mpc~\ and use 
the AB photometric system (lOke fc GunnllT98^ . 

2. THE SNLS DATASET 

This section describes the observational dataset that 
we use for our analysis in this paper. We first describe 
the dataset of SNe la, and then the datasets of both the 
SN la host galaxies, and the general galaxy population. 

2.1. Supernova data 

Our SN la data come from the Supernova Legacy Sur- 
vey (SNLS). The SNLS is a "rolling" search for dis- 
tant SNe with a primary science goal of using 700 high- 
redshift SNe la to determine the average equation-of- 
state parameter of dark energy, < w > fsee lAstier et alJ 
2006 ). SNLS exploits th e square-degree Megacam cam- 
era IjBoulade et alJl200,'^ on the Canada- France-Hawaii 
Telescope (CFHT) to conduct repeat g'r'i' z' imaging 
of 4 low Galactic extinction fields (named Dl to D4, 
see ISuHivan et all lloOfit for the field coordinates), im- 
aged as part of the "deep" component of the five year 
CFHT Legacy Survey (CFHTLS). The data are time- 
sequenced with observations conducted every 3-4 nights 
in dark time, allowing the construction of high qual- 
ity multi-color SN light-curves. Supplementary u*-band 
data is also acquired, though these data are not time- 
sequenced. SN candidates are observed spectroscopically 
at the ESO-VLT, Gemini or W. M. Keck telescopes to 
confirm their nature and to measure a spectroscopic red- 
shift eith er from host galaxy features or from the SN 
itself fsee lHowell et al.l l2005. for an overview of the spec- 
troscopic typing). A description of the real-time search 
operations and the criteria for fo llowing SN candidat es 
spectroscopically can be found in lSullivan et alJ l)2006D . 

The survey is now in its third year, and each of the 
four survey fields has now been observed for at least two 
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complete "seasons" , each spanning the 5-6 months the 
field is visible from Mauna Kea. This paper uses data 
from the first two seasons of Dl, D2 and D4, as well as 
the second and third seasons of D3. (The first season 
of D3 was performed in a pre-survey period where the 
data quality and spectroscopic completeness were sub- 
stantially lower.) During this period, 116 spectroscopi- 
cally confirmed SNe la were obtained over the redshift 
range 0.2 < z < 0.75, of which 100 were detected during 
the period of the survey for which an accurate rate ef- 
ficiency calculation can be performed (^Neill et al..i200^ 
and S I4.2|l . To this confirmed SN la sample, we add 24 
probable SNe la (see § 14 .211 which lack a spectroscopic 
identification but for which a reliable SN la photometric 
redshift can be estimated and which possess an excellent 
hght-curve fit to a SN la. This set of 124 SNe la forms 
the primary sample studied here. 

2.2. Galaxy data 

The images that we use to measure the galaxy param- 
eters are also constructed from data taken as part of the 
CFHT-LS. Deep optical stacks were generated for each 
filter (u *q'r'i'z') in each field from "Elixir" processed 
images (jMaenier &: Cuillandrc 2004) available from the 
Canadian Astronomy Data Centre (CADC)^'^. We sum- 
marise the main details of our processing steps here. A 
precise astrometric solution, accounting for distortion, is 
assigned to every image frame, and the seeing and photo- 
metric quality determined by flux measurements of ter- 
tiary standard stars in the CFHT-LS fields compiled by 
the SNLS team. Any two-dimensional sky variation is 
removed from each frame by fitting the background spa- 
tially and subtracting the resultant fit. Each individual 
frame has a weight-map associated with it, containing 
the uncertainty in each pixel from considerations of pho- 
ton noise from the sky background and object photons. 
Known bad pixels (as determined from the Elixir flat- 
fields) and saturated pixels are assigned a weight of zero. 

For the SN host galaxies, we require deep optical stacks 
with no SN light present when we measure the fluxes of 
a given object (the presence of SN light would otherwise 
bias the galaxy flux measurements). We therefore gen- 
erate stacks on a "per-season" basis; stacks for season 
one SNe are generated using only data taken in season 
two and later seasons, and so on. The period between 
observing seasons for each field - 6 to 7 months - allows 
ample time for SN light to have faded to an negligible 
level. Every image is resampled to a common pixel co- 
ordinate system, excluding images that do not meet the 
data quality criteria for a given field/filter/season com- 
bination. The goal is to maximize exposure depth whilst 
retaining excellent seeing in the final stacks (hence dif- 
ferent seeing cuts were used for each field/filter /season 
to account for varying observing conditions at different 
times of each year). We also require that for a given 
field/season combination, the stacks in each filter have a 
similar seeing. The typical seeing of the final stacks is 
0.7-0.8". The SWARF package (version 2.15.6) with a 
LANCZ0S3 kernel is used for the resampling. The re- 
sulting resampled frames are combined using a weighted 
average, with a sigma-clipping to remove artifacts (from 
satellite trails and cosmic-rays) using a custom-written 

http://cadcwww.dao.nrc.ca/ 



routine. 

We also require flux information for the general galaxy 
population to act as a comparison sample to the SN la 
host galaxies in later sections. For this we use data from 
the two fields for which the deepest data exists (Dl and 
D4), particularly in u* . We use the same stacking al- 
gorithm as above, but produce stacks including all the 
data for a given field - the field galaxies have no contam- 
inating SN light to concern us - hence these stacks are 
deeper than those used for the host galaxies. 

Photometric zeropoints are derived from observations 
of standard stars taken from ISmith et a l. (2002), and 
transformed to the Sloan Digital Sky Survey (SDSS) pho- 
tometric system; these zeropoints are determined sepa- 
rately for each of the deep fields in the month in which 
the photometric reference epoch is defined. Two of the 
fields (D2 and D3) are in the SDSS, and this provides an 
independent check of the zeropoints. The CFHT-LS ob- 
serves in a filter system close to the ISmith et all l|2002D 
u'g'r'i' z' SDSS system, though the Megacam u* filter 
differs from the SDSS u' filter, being designed to take 
advantage of the improved UV capabilities of CFHT and 
Megacam. Ou r effec tive filter responses can be found in 
ISullivan et all ll20(B . 

2.3. Host galaxy identification and flux measurement 

The correct host galaxy for every SN la in our sample is 
identified as follows. We begin by measuring the SN po- 
sition from point-spread function (PSF) fits to the SN in 
subtraction images on epochs when the signal-to-noise is 
largest. These SN positions are accurate to a fraction of a 
pixel, or ~ 0.02" (one megacam pixel is 0.186"). By con- 
verting these pixel positions to RA and DEC, we can then 
precisely identify the location of the SN on the u*g'r'i'z' 
deep stacks even though no SN light is present, as all the 
frames share a common astro metric system. For each 
SN la, we use SEx tractor 2.4.4 l|Bertin fc Arno uts 1991, 
see lHolwerdal2005l for an excellent user guide) to identify 
all potential host galaxies within 10" of the SN position, 
and simultaneously measure photometric information on 
the host galaxies. SExtractor is used in dual-image mode 
for the photometric measurements; detections are per- 
formed in the «' filter (the filter with the deepest data) 
and measurements performed in each of u* g'r'i' z' . 

Various photometric measurements and galaxy struc- 
tural parameters are recorded for each galaxy. We use 
the SExtractor MAG_AUTO flux measures and associ- 
ated errors throughout this paper. The MAG- AUTO 
apert ure is a flexible elliptical aperture (e.g. iKronI 
11980(1 with a characteristic "Kron" radius; we con- 
figure SExtractor to measure fluxes inside 2.5 Kron 
radii. Analytically, this results in '^90% of an object's 
light being measu red in each aperture l|Infantel Il987t 
iGraham fc Driver! 12005). though this may break down 
for very faint, highly concentrated galaxies at the de- 
tection limit (Bernstein et al. 2002; Benftcz et al. 2004 
see discussion in IGraham fc Driverl I2005ll . The use of 
SExtractor in dual-image mode (and the similar PSFs of 
the stacks in the different filters) ensures that the same 
size MAG_AUTO aperture is used for the different mea- 
surements of a given galaxy, minimizing aperture mis- 
matches. Statistical errors in each flux measurement are 
estimated using the weight-image of the final stack as an 
r.m.s image in SExtractor. 



4 



Sullivan et al. 



The identification of the correct host galaxy is not al- 
ways straightforward - occasionally the closest host in 
arcseconds is probably not the correct identification. In- 
stead we use a separation normalized by the apparent size 
of the galaxy to which the SN is being compared (Fig- 
ure^. We calculate the separation of the SN from each 
candidate host galaxy in terms of the elliptical radius (i?) 
along a line connecting the SN pixel position (xsn,ysn) to 
a given host center (a;gai,j/gai)- The elliptical shape is de- 
termined by SExtractor, defined by semi-major (ta) and 
semi-minor (r^) axes together with a position angle (6*), 
with R given by 

— CxxX^ + CyylJ^ -|- CxyXrUn (1) 
where Xr = Xsn-Xgal, Hr J/sn-t/gal, C'xx = C0s'^i9)/r\ + 

sin^(6')/r|, Cyy = sin^(6')/r^ + cos'^ (0) / r^, , and Cxy = 
2 cos{9) sm{9){l/r\ + ^/r%). R is determined for every 
candidate host galaxy, the hosts ordered by R, and each 
SN is assigned to the host that is nearest in terms of 
this parameter. The isophotal limit of a given object 
corresponds to i? ^ 3; we consider only galaxies with 
R<5. In some cases (~7%), no host galaxy is identified 
by SExtractor with i? < 5; we then measure the flux and 
flux error inside a 3" diameter aperture centered on the 
SN position, and use this as the "galaxy" flux in what 
follows. Figure n shows a visualization of the technique 
for four SNe la. 

We also measure the properties of the general galaxy 
population using the deeper stacks constructed as de- 
scribed above. We use SExtractor with the same param- 
eters as for the SN host galaxies. 

3. GALAXY PROPERTIES 

This section details the techniques we use for convert- 
ing the galaxy properties into characteristics that can be 
used for our subsequent analysis in § |S1 and § El We first 
discuss the technique of converting the observed fluxes of 
S I2.3I into galaxy stellar masses and star-formation rates, 
both for the host galaxies themselves and for the general 
field galaxy population. We then detail a number of con- 
sistency checks we perform on these derived parameters. 

3.1. Galaxy SED fitting 

In order to investigate SNe la in relation to the mean 
properties of the stellar population from which they 
were formed, we need to derive properties such as stel- 
lar mass and star-formation rate (SFR) for the galax- 
ies. To do this we fit a series of template galaxy 
spectral energy distribution (SEDs) to the broadband 
fiuxes available for each galaxy (§ I2.3|l . and then use 
the best-fit SED for the estimation of the variou s 
properties (see for example iBrinchmann fc EUid l200(]j) . 
This technique is similar to that employed by pho- 
tomctric r edshift cod e s (e.g . 'Gwvn & HartwicW 119961; 
[Bolzonella et al.l 120001 : iLe B orgnc & Rocca-Volnicrange 
|2002|) . The best-fit SED is determined using a mini- 
mization between the observed fluxes, the corresponding 
flux errors, and the synthetic photometry generated by 
integrating the template SEDs through the SNLS effec- 
tive filter responses. For the host galaxies, the spectro- 
scopic redshift is known from the SN confirmation spec- 
trum and held fixed, reducing the uncertainty in the de- 
rived properties of each galaxy. For the field population 



no spectroscopic information is available, and the red- 
shift is left as a free parameter in the fits. 

Our set of synthetic templates is computed 
with the PE GASE.2 galaxy spectral evolu- 
tion^ llFioc & R occa-Volnicrangc 199'3 

Le Borgne fc Rocca-Volmerange ,,2002: .Le Borgne et al] 
2004|) . This code and the SEDs that it gener- 
ates have been used extensively in the literature 
to co nstrain the prop erties of high-reds hift galaxies 
fe.g. iGlaze brook et alJ I200I iMcCarthv et'al.1 I200I 
iGrazian et al.i i200ffl ! We use eight scenarios which 
evolve self-consistently with age; details of the physical 
parameters defining them can be found in Table 1 of 
Le Borgnc & Rocca-Volmerangc (200^. We assume a 
Krouoa (2001) initial mass function (IMF), and assume 
the IMF is unive rsal across environment (see discussion 
in lKroupall2002|) . The synthetic SEDs are the sum of 
light emitted by stars and nebular emission (continuum 
and lines), including attenuation by dust with a King 
model or a plane-parallel slab geometry. The quantity of 
dust evolves consistently with the amount of gas present 
in the galaxy. 

An SED is computed at 69 time-steps in each of the 
8 scenarios, giving a total of 552 template SEDs. When 
fitting a given galaxy, we use only templates younger 
than the age of the Universe at the redshift of the galaxy. 

3.2. Galaxy derived properties 

Two physical parameters, for both the SN la host 
galaxies and for the general galaxy population, are of par- 
ticular interest in this study: the galaxy mass, and the 
galaxy SFR. The first is the total stellar mass of a galaxy, 
the total current mass in stars of all types and ages. This 
is derived by integrating the total star-formation history 
(SFH) of the best-fit scenario up to the best-fit age, and 
subtracting the mass of stars that have died. 

The second parameter is the amount of recent star- 
formation that the galaxy has experienced. There are 
several approaches that could be taken. First, one could 
simply divide the ultra-violet (UV) luminosity of each 
galaxy by a conversion factor to obtain a UV-derived 
SFR fe.g. Eonas et aUl987|) . However, such an approach 
tends to over-estimate the SFR in old stellar systems due 
to the contaminating presence of evolved stars. A second 
approach to determining instantaneous SFRs is to use 
nebular emission lines (Kennicutt 1998). Unfortunately, 
only very limited information on the appropriate lines is 
available for our host sample from the SN confirmation 
spectra, with large uncertainties from aperture losses, SN 
light contamination and flux calibration. 

We instead estimate the mean SFR from our best- 
fitting scenarios, averaging the SFR over a longer time 
interval. This provides an automatic correction for the 
UV light from old evolved stars. To define the period over 
which the star-formation should be averaged, we carried 
out simulations where the mean SFR was estimated from 
synthetic u*g'r'i' z' photometry measured both on the 
idealized models and on models with stochastic SFHs 
(see S I3.3|l . These simulations showed that mean SFRs 
on periods of 0.5 Gyr can be reliably recovered without 
significant systematic errors. On shorter timescales sys- 
tematics can be introduced, particularly for galaxies for 
which the redshift is not known (see Tableland discus- 
sion in ij l3.3|l . 
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For the general galaxy population, an alternative ap- 
proach would be to use a parameterized form of the mass 
or luminosity functions, either directly in ferred from 
high- r edshift spectroscopi c surveys ^e.g.. |Ilbert et all 
1^05^ 'Wi Umer et aL * "200^: ' Bundv et all *2005') or low- 
redshift functions adjusted for evolution to high redshift. 
While this is possible for the global mass function, no 
equivalent distribution exists for the SFR property of 
field galaxies, requiring a conversion from (for example) 
a rest-frame [/-band or UV luminosity function. Fur- 
thermore, since a goal of this work is to investigate host 
galaxy properties binned by the type of the host galaxy 
(§Ellj mass and luminosity functions derived from galax- 
ies categorized in the same way would be required, and 
these are not currently available. 

Uncertainties in the derived parameters for our galax- 
ies arise from both statistical and systematic errors (dis- 
cussed in the next section). Statistical errors derive from 
the photometric measurements of the galaxies, and are 
accounted for by considering the full range in the qual- 
ity of the 552 template fits (as defined by the statis- 
tic). This maps out a probability surface in the stellar 
mass and SFR parameters determined from the fits, and 
allows estimates of the statistical uncertainties of each 
of the derived physical parameters; lower signal-to-noise 
photometry will be consistent with a wider range of fit 
templates and hence a wider range of stellar masses and 
SFRs. Invariably, the derived masses and SFRs possess 
non-symmetrical error-bars, particularly for the fainter 
systems. In some cases, fainter galaxies can have more 
than one minimum in the space. In this case, the 
best-fit is considered as the solution, but the error bars 
in mass and SFR are extended to the largest interval of 
uncertainty, covering the range of parameters between 
the minima. The errors are therefore conservative for 
galaxies with several minima. 

3.3. Systematic errors 

Systematic errors in our derived galaxy stellar masses 
and mean SFRs potentially arise from many sources. 
These include the wavelength range of our input photom- 
etry, and our choice of spectral libraries, IMFs and SFHs 
that form the template SEDs used for the fitting, and are 
of course harder to estimate. We first test the accuracy of 
the photometric redshift estimates by running fits to the 
host galaxies (with known spectroscopic redshifts) with 
the redshift left as a free parameter. Figure|2](left) shows 
these spectroscopic and photometric redshift estimates 
for the SN la host galaxies. For galaxies with SExtrac- 
tor detections, the median of |Az| = | ^spec — -Zphot I is 0.02 
(the 90%-ile is 0.15), and the median of |Az|/(l -I- Zgpoc) 
is 0.012 (90%-ile is 0.09). These measures contain pho- 
tometric redshifts from galaxies with a significant pho- 
tometric uncertainty due to their faintness (and hence 
uncertainty in the photometric redshifts). To confirm 
that these error estimates are reasonable, we calculate 
the dispersion |Az|/cr, where cr is the appropriate er- 
ror in the photometric redshift. 90% of the photometric 
redshift estimates lie within 2(j of their corresponding 
spectroscopic redshift. 

We next compare the host galaxy properties derived 
from the template fit with a known redshift against those 
de rived from the pho tometric redshift fit (c.f. figure 3 
of IBundv et all 12001 . Such a comparison will show 



whether the photometric-redshift properties of the galax- 
ies we study can reliably trace the true galaxy properties 
in the mean. We compare the stellar mass and mean 
SFR estimates for galaxies using the spectroscopic red- 
shift with the same estimates when the redshift is allowed 
to float (Figure OJ. The agreement for both quantities is 
reasonable, with the mass estimates appearing more ro- 
bust. No mass or SFR-dependent trends are seen. The 
median difference AM = Mgpoc-Afphot is 0.026 dex (90% 
of AAf lie within 0.2 dex), and the difference ASFR is 
0.033 dex (90% of ASFR he within 0.38 dex). A compar- 
ison of the mass function derived from our photometric 
redshift fits to published mass functions of field galaxies 
will be presented in an upcoming paper (Le Borgne et 
al., in prep). 

The next potential source of concern is the lack of 
near-IR data and the impact on our estimates of stel- 
lar masses, and the lack of very short wavelength data 
and the impact on the derived SFRs. These could be 
particularly pronounced in systems with both a young 
and an old stellar populations, such as those galaxies 
experiencing recent star- formation events. Clearly the 
uncertainties in derived stellar masses and mean SFRs 
will be larger in the absence of UV and near-IR data; 
however as long as the error-bars in these properties re- 
flect this then this additional uncertainty will be carried 
through in our analysis. A larger concern would be any 
systematic under- or over-estimation of the stellar masses 
and SFRs. 

We investigate this by using our standard templates 
to fit photometry generated from SEDs with stochastic 
SFHs. Synthetic SEDs are formed from the composite of 
3 random exponential SFHs each with a different age [t) , 
mass fraction and r. One SFH is constrained to be old 
(t>lGyr), one of intermediate age (200Myr<i<2Gyr) 
and one young (30Myr<i<100Myr). Each of these 3 
random SFHs is converted into an SED at the randomly- 
selected t, and the 3 SEDs added together to form one 
composite SED. We place this composite SED at a ran- 
dom redshift (0.2<2:<0.75), generate synthetic u*g'r'i'z' 
photometry, and see how well PEGASE.2 recovers the 
stellar mass and mean SFR using the 8 idealized scenar- 
ios of § 13.11 both when the redshift is fixed and when it 
is left as a free parameter. 

The results are given in Table ^E^nd are encouraging. 
For the stellar masses, we see only small differences in the 
mean between the input and recovered value with a stan- 
dard deviation of around 0.3dex. For the mean SFRs, 
again there is no significant systematic offset though the 
standard deviation of the differences is larger particularly 
when the redshift of the galaxy is not known. 

We also briefly investigate the effect of dust on our stel- 
lar mass and mean SFR estimates. Though dust is in- 
cl uded in the idealized scenarios (see 13. H and discussion 
in lLe Borgne fc R,occa-Volmerang'3l2002yi . this will only 
provide an average correction to our derived properties. 
We experimented with adding extra dus t to our simu- 
lated SEDs (up to E(B-V)=0.5 using a iCalzetti 6^11 
120001 extinction law). As expected, this can lead to a 
systematic under-estimation of the mean SFR (up to 0.4 
in dex in some cases) and a small increase in the stan- 
dard deviation, though the accuracy of the stellar mass 
estimates was not affected. 
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Finally, we estimate systematic uncertainties by de- 
riving the galaxy properties using a similar technique 
but different templates. We use the photometric redshift 
code of [Gwvn & Hartwick (1996), which uses a differ- 
ent set of (empirical) gal axy templates ( Colema n et alJ 
ITosnt iKinnev et alUlOOl to those used in Z-PEG. We 
determine masses by fitting iBuzzonil ()2005j) population 
models to the best-fitting templates which results in a 
mass-to-light ratio which can be used to convert the 
galaxy luminosities into masses. SFRs are determined 
by divi ding the f lux at 2800A by 4.8 x 10^^ erg s~^ Hz~ ^ 
ijBuzzo ni 200^ iMadau. Pozzetti. k DickinsonI Il998hj) . 
with a small correction for early-type SEDs to account 
for the contribution from older stars (see Buzzoni 2002) 

The two different techniques for estimating the masses 
and the SFRs showed a good agreement. The mass and 
SFR estimates agree to better than a factor of two - 
mean difference in mass was 0.35 dex (in the sense Z- 
PEG measured smaller masses) with an r.m.s. scatter 
of 0.24, and the mean offset for the SFR estimates was 
0.15 dex (in the sense Z-PEG measured larger SFRs) 
with an r.m.s. scatter of 0.29. Most importantly, no 
mass-dependent or SFR-dependent trends were detected. 
Given that Z-PEG is used to calculate masses and SFRs 
for both field and host galaxies, these offsets in mass or 
SFR will cancel in our analyses. 

4. INCOMPLETENESS CORRECTIONS 

In the next two sub-sections we discuss incompleteness 
in our sample. We first discuss the incompleteness of our 
field galaxy population due to the limiting magnitude 
of our deep field stacks, an effect common to all galaxy 
redshift surveys. We then discuss the incompleteness of 
our SN la sample. 

4.1. Galaxy incompleteness corrections 

Incompleteness will affect our stellar mass and mean 
SFR distribution functions of the field galaxy popula- 
tion. Our survey is magnitude limited; galaxies with a 
given absolute magnitude (and spectral type) will be- 
come fainter than this magnitude limit at different red- 
shifts, which may be less than the largest redshift we 
consider here (z = 0.75). Wc adopt the traditional Vmax 
method (e.g. iSchinidt 1968: Fcltcn 1976,) to correct for 
this effect. Vmax is defined as the co- moving volume 
within which each galaxy - as defined by its absolute 
magnitude and fc-correction from best-fitting SED tem- 
plate - would remain in our sample i.e. within which 
it would satisfy the limits in apparent magnitude of the 
current depths of the optical stacks. Each galaxy is then 
weighted by a factor T^urvcy/Knax when computing the 
various mass and SFR distributions used in § |5l where 
Ksurvcy is the total survey volume. 

The SExtractor photometric measurements are per- 
formed by detecting in the i' filter, and the depth in 
this filter defines the limiting magnitudes. The limiting 
magnitudes on the i' Dl and D4 stacks are determined 
by inserting fake sources into the stacks, and measur- 
ing the fraction recovered by SExtractor as a function of 
magnitude; the limiting magnitude is that at which 50% 
of fake sources are recovered. 

4.2. SN incompleteness and rate calculation 



Incompleteness will affect our SN la sample in several 
ways. Firstly, we could miss SNe la altogether due to 
inefficiencies in our search pipeline (the detection effi- 
ciency). Secondly, detected SNe la in the redshift range 
of interest could be lacking a spectroscopic observation 
because such an observation would be too challenging, 
for example due to host galaxy contamination or the SN 
maximum light falling during a bright-moon period. Fi- 
nally, SNe la could lack a spectroscopic confirmation be- 
cause poor weather affected the spectroscopic scheduling. 

To corre c t for these effects, we adopt the scheme of 
iNeill et al.l l|2n0fit) . who calculated the SNLS SN la rate 
over 0.2 < z < 0.6, but we extend the calculation to z = 
0.75. The first source of incompleteness, detectability, 
is accounted for using simulated SN data. Fake SNe are 
placed in real SNLS data, and the recovery efhciency as a 
function of SN brightness, position, sky background, host 
brightness, expos ure tim e, seeing, and sky transparency 
is calculated (see iNeilFet al. 2006). 

We also determine the observing window within which 
maximum-light of a SN la would have to occur in order 
to be considered for a spectroscopic observation and for 
which we can derive a reliable photometric classification 
(see below). SN la candidates detected at the start or 
end of an observing season lack a full light-curve and 
therefore are not usually observed spectroscopically as 
any resulting light-curve fit would be far more uncertain. 
For the purposes of this study, these cut-offs are defined 
as follows, i) The SN must have at least 2 i' and at least 
1 r' observation between -15 and -1.5 days in the SN rest- 
frame, ii) There must be at least 1 g' observation between 
-15 and -1-5 days, and iii) There must be at least 1 r' or i' 
observation after -1-11.5 days but before -1-35 days. These 
criteria cull SNe la whose light-curves do not properly 
sample maximum light and for which a light-curve width 
measurement would be correspondingly more uncertain, 
and for which a color near maximum was not measured. 
Any confirmed SNe la which do not meet these criteria 
are excluded from this analysis; this reduces our sample 
from 116 SNe la to 100. Our simulations then provide, 
for a given field observing season, the efficiency required 
to convert an observed "per season" rate (the number of 
SNe la with 0.2<z<0.75 meeting the light-curve criteria 
defined above) into a yearly rate. These efficiencies can 
be found in Table H 

The spectroscopic incompleteness - the fraction of can- 
didate SNe la detected but never observed spectroscop- 
ically - is the most challenging source of uncertainty to 
address in this study. We assess the incompleteness us- 
ing the photometric SN selection method presented in 
l^uUivan ct al. (2006). This technique fits SN hght -curves 
in the absence of a spectroscopic redshift, and returns the 
best-fitting SN la parameters (redshift, stretch, E^^^y, 
and a dispersion in the SN peak magnitude, dm), as well 
as a guide that the candidate under study is a SN la. 
The software was run on g'r'i'z' photometry for all SN 
candidates (including those observed and typed spectro- 
scopically) discovered during the period covered by the 
simulations above. A comparison of SN spectroscopic 
redshift and SN photometric redshift for all SNLS SNe la 
can be found in Figure|21 (Note that as the SN photomet- 
ric redshift code makes an assumption about the cosmol- 
ogy when fitting the light-curve, the resulting photomet- 
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ric redshift cannot themselves be used for determination 
of the cosmological parameters.) The mean difference 
^spec — -Zphot for ah SNe la is -0.007; the standard devi- 
ation is 0.088. For the SNe la included in this study - 
over 0.2<z<0.75 and which meet the light-curve criteria 
- the standard deviation is 0.077. Occasionally, a spec- 
troscopic observation of a SN candidate yielded a spec- 
troscopic redshift but no definitive type; in these cases 
we still allow the redshift to float for comparison with 
the spectroscopic redshift. 

There were 286 SN candidates without a spectroscopic 
confirmation with a SN la photometric redshift in the 
range of interest. 80 do not meet the light-curve coverage 
criteria that we apply to the spectroscopically confirmed 
sample, and are excluded. This ensures that all candi- 
dates have a reliable stretch and g' observations. We also 
exclude 2 SN candidates where the spectroscopic red- 
shift was in disagreement with the photometric redshift 
where Zphot > ^spcc + 0.135 x (1 + Zspcc) (FigureEJ. In 
these cases, the SN spectra were ambiguous and no type 
could be determined; usually this was because a SN Ib/c 
spectrum provided a similar quality fit to the observed 
spectrum as did a SN la. 

The remainder (204) were then culled based on their 
light-curve parameterization and the of the light-curve 
fits in the different observed filters. We exclude 70 SN 
candidates with a fitted stretch >1.35 (very effective at 
removing SNe IIP), and 75 candidates whose per de- 
gree of freedom in the fit was >10. None of our spectro- 
scopically confirmed SNe la had fit parameters in either 
of these ranges (see Figure QJ. 

These simple cuts leave 67 candidate SNe la. A further 
visual inspection of the remaining candidates revealed 
that although many of these had acceptable overall 
fits, they were too blue in g' — r' before -1-5 d when com- 
pared to a SN la template. Hence, our final statistical 
cut removes 35 objects with a poor g' fit at early-times 
(-15<d<-|-5), or with a g' mean dispersion 

\ ^ yyobs -9 m 



where Ng is the number of g' light-curve fiuxes gobs with 
error g^rr over -15<d<-|-5, and gmodei is the SN la tem- 
plate model. This cut does not remove SNe that are too 
red (negative dispersions), which could be indicative of 
extinction, only those that are too blue. The various 
culls that we use are summarized in Table |21 This leaves 
24 SN la candidates to add to the 100 spectroscopically 
confirmed SNe la sample; details can be found in TableEl 
For these SNe, the redshift of the SN la is taken to be the 
SN la photometric redshift where no spectroscopic red- 
shift was available. The impact on our results of these 
photometrically classified SNe la is discussed in § El 

5. THE SUPERNOVA lA RATE AS A FUNCTION OF 
HOST GALAXY PROPERTIES 

In this section, we examine the SN la rate as a func- 
tion of the properties of their host galaxies, and use these 
rates to place constraints on the composition of the to- 
tal SN la rate. We first classify our SNe la into three 
sub-groups based upon the nature of star-formation in 
their host galaxy. The first group comprises SNe la lo- 
cated in host galaxies with a zero mean SFR from the 



SED fits ("passive" galaxies). We then use the specific 
star-formation rate, sSFR, defined as the star-formation 
rate per unit stellar m ass fe. g. iGuzman et aD Il997t 
iBrinchmann fc EllisI 1200 0: Brinc hmann et alJl2004|) . to 
classify the star-forming SN la host galaxies. With units 
of yr^^, sSFR is essentially a measure of the inverse of 
the formation time-scale for a given galaxy: high sSFR 
galaxies will form the mass in their stellar populations on 
shorter times than low sSFR galaxies. The second group 
of hosts, defined as -12.0 < log(sSFR) < -9.5, have 
a small or moderate amount of star-formation relative 
to their stellar mass, and therefore are likely to possess 
a substantial evolved stellar component as well as young 
stars. The third and final group, with log(sSFR) > —9.5, 
have a large amount of star-formation relative to their 
stellar mass, and therefore stellar populations that have 
a large component made up of young stars. Broadly 
speaking, the second group tends to comprise normal 
star-forming galaxies such as the Milky Way, whereas 
the third group tends to include vigorously star-forming 
and lower mass dwarf galaxies. The division is illustrated 
in Figure[5l which shows the distribution of all the SN la 
host galaxies in the stellar-mass/ SFR plane. 

An attraction of the SED-fitting technique is that 
galaxies can be classified according to their star- 
formation properties without regard to morphology. 
This has some advantages; in morphologically-selected 
samples of distant galaxies, a significant fraction of 
spheroidals have been shown to possess both "blue 
cores" and weak O ii emission lines in their spectra 
Ipylenantcau et al. 2001; Treu et al. 2005), interpreted as 
a signature of young stars in these galaxies. There is 
also some evidence for recent star-formation in some 
early-typ e galaxies from the near-UV color-magnitude 
relation ijYi et alJ l2005V The advantage of the effective 
color-selection of SED-fitting is that rather than assume 
all early-type galaxies consist purely of old stars, the 
evolutionary model fitting places no a priori constraint 
on the type of SFH that can be fit to a given galaxy. 
Young populations are perfectly possible in morpholog- 
ically spheroidal galaxies, and hence any young stellar 
population which has a significant impact on a galaxy's 
colors will be reflected in the best-fitting galaxy SFH. 

We examine how the SN la rate varies as a function 
of this specific SFR in Figure El The field galaxies are 
binned according to their value of specific SFR, and the 
total field-galaxy stellar mass in each bin of specific SFR 
is calculated. We use the Knax technique of S 14.11 to in- 
completeness correct this distribution. We then bin the 
SN la host galaxies by specific SFR in the same way and 
calculate the number of SNe la per unit stellar mass as a 
function of the specific SFR of their host galaxies. In Fig- 
urejni the rate in passive galaxies, where the SFR is zero, 
is shown as a hashed area starting at log(sSFR) — —12 
(the height of this hashed area represents the statistical 
uncertainty in the measurement). 

An increase in the rate of SNe la per unit stellar mass 
with increasing specific SFR of the host galaxy is clear. 
The difference between the rate in passive galaxies and 
the most vigorous star-forming systems is about a fac- 
tor of 10; furthermore, the increase in the rate is a 
fairly smooth function of the specific SFR. The general 
trend of Figure El can be compared to that observed by 
iManmicci et alJ l|2005D in the local universe using a mor- 
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phologically classified sample of local SNe la host galax- 
ies. We illustrate this in Figure El by over-plotting the 
Mannucci et al. data on the SNLS results; the Mannucci 
et al. evolution is shown normalized to the SNLS rate in 
passive galaxies. The trend is very similar, though there 
is the obvious caveat that the link between specific SFR 
and galaxy morphology is not a straightforward one-to- 
one mapping. 

As noted by other authors, this relationship is difhcult 
to reconcile with a model for SNe la that originate s solely 
from a n old evolved stellar population. Man nucci et al.l 
l|20n(it) and SB05 instead model the SN la rate (SNRia) as 
a composition of two separate components: a "prompt" 
component, with a short delay time, and an "old" com- 
ponent, with a long delay time. The most general form 
for the SNRia as a function of time is simply the con- 
volution of the SFR, Mnew , and the probability function 
for getting a SN la from a stellar population of age t, P, 
i.e. 

SNRia(0=/ Mnen,{t')P{t-t')dt' (3) 

Jo 

This rate can be simply modeled by making the assump- 
tion that P can be well represented by two components. 
One has a peak of B at time i = and is zero at all other 
times (this represents a very short delay time) , the other 
has P — A constant with time (and represents long delay 
times), i.e. 

SNRia(t)=A/ Mnew{t)dt + BMnew{t). (4) 

Jo 

As Mtot(i) = Mnewit)dt, where Mtot(0 is the total 
mass of a galaxy at time t, this equation models the 
probability of a SN la exploding in a given galaxy as 
depending on both the mass and the instantaneous SFR 
of that galaxy. A and B are constants which relate the 
total mass and the SFR of a galaxy to the SNRia in 
that galaxy (which SB05 fix using the observations of 
[Mannucci ct al. (2005)). In effect, A is the SNRia per 
unit mass of the old component, and B the SNRia per 
unit SFR of the young component. The model predicts 
that SNRia is linear in both host galaxy mass and SFR. 

Note that the definition of mass above is slightly dif- 
ferent to that measured for the SN host galaxies § 13.11 
For the hosts, we measure the total stellar mass i.e. the 
total mass currently in stars in each galaxy. Mtot above 
is the integral of the SFH for each galaxy ~ no correc- 
tion is made for stars that have lost mass at the end 
of their stellar evolution. The numerical differences in 
these mass definitions are shown in Figure [7| with the 
differences being largest in older, lower SFR systems. 

The model of equation is a simplification of the real 
physics. For the prompt component, the model implies a 
zero delay-time between star-formation and SN la explo- 
sion; some non-zero delay-time to account for main se- 
quence lifetime and subsequent accretion onto the white 
dwarf is obviously required. For the old component, the 
equation simplifies the complex relationship between the 
SN la delay-time and the age of the stellar population by 
using a simple constant probability rather than a more 
complex exponential or Gaussian delay time distribution. 
In reality, for a coeval population, after a few billion years 



the probability of a SN la will likely decrease as the stellar 
population ages and fewer progenitor stars are available; 
this could cause an over-estimation of the SN la rates in 
the oldest, passive stellar systems. 

Yet, the se approximations may not be that poor. 
iMannuccTet al. (2006) find that the local SN la delay- 
time distribution is well represented by a prompt compo- 
nent modeled as a Gaussian centered at i = 50 Myr plus 
a component modeled as an exponential with a decay 
time of 3 Gyr. These two terms can be approximated by 
a delta function at time t — plus a constant probabil- 
ity thereafter. Furthermore, this equation parameterizes 
the rate in a convenient form relating to galaxy proper- 
ties that are relatively straightforward to measure using 
the SNLS dataset introduced above. 

In the next sections, we illustrate this model by con- 
sidering the SN la rate as a function of both host galaxy 
stellar mass and host galaxy SFR taken separately, and 
then constrain the A and B parameters using a bivariate 
fit to the mass/SFR data. 

5.1. SN la rate as a function of host galaxy stellar 

mass 

We first attempt to separate any component of the 
SN la rate which may depend on the stellar mass of 
a galaxy, from any component which depends on the 
amount of recent star-formation i.e. to separate the two 
components of equation Q. We can do this using our 
sample of SNe la which exploded in passive galaxies, and 
which in our models have zero recent star-formation ac- 
tivity. By binning the SN la host galaxies according to 
their mass, and dividing by the equivalent (incomplete- 
ness corrected) distribution of the general galaxy popu- 
lation in each mass bin, we can examine how the proba- 
bility of a SN la explosion is related to the mass of the 
host galaxy. 

We perform this comparison in Figure |H1 where we 
show the dependence of the SN la rate as a function 
of mass of the host galaxy. To account for the (invari- 
ably non-symmetrical) errors in the mass determinations 
of the SN host galaxies, we have performed a Monte- 
Carlo simulation with 5000 realizations of each SN host 
with the masses for each host drawn from the estimated 
probability distribution for that host. We then bin this 
Monte-Carlo population, normalize to the total number 
of hosts in the sample, and use this binned distribution 
in our analysis. Such a procedure produces a distribution 
which accounts for observational statistical uncertainties, 
which can be considerable in the fainter host sub-sample 
where some of the measurements effectively provide only 
limits on the host parameters. 

One feature of Figure|Hlis the increase in the SN la rate 
with increasing host galaxy mass, present among host 
galaxies of all types. As the SFR of the passive galaxies 
is zero (this defines their selection) , the contribution from 
the B term in equation Q is also zero in these galaxies. 
Hence the best-fitting line in log-log space should have a 
slope of 1, if equation Q is a good approximation - the 
model is linear in mass. For passive galaxies, we find that 
the slope, rimass, is fimass — 1.10 ±0.12 with a reduced 
of ~ 1.14 (if only spectroscopically confirmed SNe la are 
used in the fits, n^ass = 1-04 ± 0.13). This implies that 
the relationship between the SN la rate and galaxy mass 
is consistent with being linear and the model provides 
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and adequate fit to the data. Using this technique we 
estimate A = 5.l± 1.2 x lO^^'' SNe yr"! Mq^ . 

The relationship in later-type, star-forming galaxies is 
different. The best-fitting slopes are statistically consis- 
tent: 0.66 ±0.08 (reduced = 0.77) and 0.74 ±0.08 (re- 
duced = 0.94) for the low and high specific SFR galax- 
ies respectively. This translates to an excess of SNe la in 
low mass star-forming galaxies compared to the passive 
galaxies, but a similar number in each at the most mas- 
sive end. This is as expected if the model of eqn. is 
correct. Doubling the mass in passive systems will dou- 
ble the SN la rate, but doubling the mass in star-forming 
systems will only double the rate if the SFR doubles as 
well. As the most massive star-forming systems generally 
have lower specific SFRs than the lowest mass systems, 
the increase in the SN la rate with host mass in star- 
forming systems is not linear. The result is that while 
the prompt component SNe dominate in low-mass star- 
forming galaxies, the old component is more important 
in higher-mass star-forming galaxies. 

5.2. SN la rate as a function of host galaxy mean SFR 

Having examined the SN la rate in galaxies dominated 
by old stellar populations, we now examine the SN la rate 
in very young populations via the SFR, which we deter- 
mine as the mean SFR over the last 0.5 Gyr of a galaxy's 
SFH as determined from the best-fitting template sce- 
nario fS I3.1|l . We first bin both the SN la hosts and the 
general galaxy population by their mean SFR. We in- 
completeness correct the general galaxy population, and 
in each bin of mean SFR for this general galaxy popu- 
lation, we sum the total stellar mass in that bin, and, 
using a value of A determined from the passive galaxies 
in § 15.11 calculate the number of SNe la expected from 
any old component. This old component is subtracted 
from each bin, leaving an excess of SNe la above that 
predicted from the stellar mass of the galaxies. We show 
this distribution in FigureEl In this figure, we have com- 
bined the SNe la occurring in all types of star-forming 
galaxies rather than the two separate star-forming pop- 
ulations shown in Figure |H1 As in § 15.11 a Monte-Carlo 
simulation is used to account for the uncertainties in the 
stellar mass and mean SFR estimates for the host galaxy 
population. 

Figure |^ shows clear evidence for a component of 
SNe la beyond that dependent on the integrated stel- 
lar mass. The number of SNe la expected in each SFR 
bin from the stellar mass is smaller than actually ob- 
served; the majority of SNe la in star-forming galaxies 
appear to arise from more recently formed stars than 
from old stars. Furthermore, the fraction of SNe la 
from the "A" component appears fairly constant with 
the SFR of the galaxy, with a range of 14% to 21%. We 
also test whether the dependence on the mean SFR of 
a galaxy is linear by fitting the slope of the line, hsfr- 
The best fit is tt-sfr = 0.84 ± 0.06 with a = 1-31 
(the slope is nspR = 0.81 ± 0.08 for the spectroscopi- 
cally confirmed SNe la). For comparison, a line with a 
slope of unity is also shown (x^ = 2.34). We estimate 
B = 4.1 ± 0.7 X 10-"* SNe yr'^ (Mq yr"!)"! using this 
approach and enforcing a linear relationship. 

5.3. Bivariate fits 



Though the last two sections provide a useful visual- 
ization of the SN la rate as a function of galaxy stel- 
lar mass and mean SFR, a more sophisticated bivariate 
technique which fits A and B simultaneously across all 
galaxy types will result in more accurate fit values, as 
the stellar masses and SFRs are partially correlated (Fig- 
ure[SJ|. We assume the relationship is a linear function of 
only A and B (i.e. nmass and nsFR are fixed), and per- 
form a generalized linear least-squares fit in the galaxy 
mass/SFR plane. This effectively fits the two compo- 
nents of the SN la rate simultaneously to the data. We 
perform this fit by converting Figure El the distribution 
of SN la host galaxies in the mass/SFR plane, into a 
probability of a SN la explosion in a given galaxy as 
a function of mass and SFR. The probability is calcu- 
lated by binning the host distribution by galaxy mass 
and galaxy SFR, and dividing the number of SN hosts 
in each bin by the (incompleteness corrected) number of 
field galaxies similarly binned. This results in a two di- 
mensional mass/SFR probability space which gives the 
likelihood of a SN la explosion in a galaxy as a function 
of both the galaxy stellar mass and mean SFR. 

We fit two forms of equation Q to this probability 
distribution: equation Q with A and B both free, and 
the same function but with B — i.e. assuming the SN la 
rate depends only on host galaxy stellar mass (clearly 
the third alternative, fixing A = 0, makes no sense as 
plenty of SNe la explode in non-star- forming galaxies). 
The addition of the B term as a free parameter (B ^ 0) 
reduces the of the fit from 67 to 37 (60 degrees of 
freedom). Performing an F-test, we find that the null 
hypothesis that the extra B term is not needed is rejected 
at >99.99% probability. The fit results when fitting for 
both A and B with nmass = ?^sfr = 1 are A = 5.3 ± 
1.1x10-1'* (i/o/70)2 SNe yr-i and B = 3.9±0.7 x 
10-^ (i?o/70)2 SNe yr-i {Mq yr"*)"! (these errors are 
statistical only). Using the formal covariance matrix of 
the fit, we calculate a correlation coefficient between the 
input masses and SFRs of -0.22, indicating only a small 
correlation. 

One important feature of Figure |51 is the slight non- 
linear relationship between SNRia and SFR. Using the 
bivariate fit method, we can also derive best-fit values 
for Timass and nsFR by examining the variation of 
the best-fit as rimass and ngpR are varied in the fitting 
function. We find minimum x^ at nmass = I-OO^qJ^q 

"■SFR = 0.98to'ii- These are consistent though more 
reliable than the estimates of § 15.11 and § 15.21 and are 
completely consistent with rimass = ?^sfr — 1- 

Clearly, in an analysis such as these systematic errors 
are likely to be at least as large as the statistical ones 
quoted above. For exa m ple, using different I MFs (e.g. 
those of lSalpeteilll955l or lRana fc Basu"199?l can vary 
the fit values of A and B b y ar ound 10-20% of their 
value compared to the IMF of lKroupal (|200lO; using dif- 
ferent spectral libraries is likely to produce changes of 
a similar magnitude. Furthermore, although dust is in- 
cluded in the PEGASE.2 models, no dust extinction cor- 
rection is made to individual galaxies, and systematic 
errors could clearly therefore be present in the SFR esti- 
mates, especially when considering the existing evidence 
for SFR-dependent extinct ion corrections iHopkins et alJ 
l200lHSullivan et al.ll200H) . Hence, the values of A and 
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B may change as more multi-wavelength data becomes 
available for the host galaxies and better constrains the 
extinction properties of the galaxies. 

5.4. Comparison to other results 

The fits of § 15.31 show that the SN la rate in galaxies 
has a linear dependence on both stellar mass and younger 
stellar populations, which we parameterize via the galaxy 
SFR. We also derived parameters A and B which allow us 
to relate the mass and SFR of a galaxy to the probability 
of it hosting a SN la. We can compare our estimates of 
A and B with other similar numbers published in the 
litera ture. 

■Manmicci et al.l l)2005f ) calculated the low-redshift 
SN la rate per unit mass as a function o f host galaxy mor- 
pholog ical type from the SN sample of iCappellaro et aTl 
l)1999tl . In low-redshift E/SO galaxies, they find a rate 
of 3.83;i;2 X 10-" SNe yr"! Mg^ (in our cosmology), 
which agrees with our determination {A — 5.3 ± 1.1 x 
10^^**) within the error-bars. Given the different tech- 
niques for determining galaxy types (SED fitting ver- 
sus morphological typing), and the different search and 
efficiency calculations, this is an encouraging level of 
agreement. Our determination of B (3.9 ± 0.7 x 10-** 
SNe yr"-'^ (Mgyr"^)"^ ) is however discrepant with val- 
ues of SB05, who derive two different values, IOI5 ^ 10"'* 
and 23tJo x 10"" in our units. If we take the SNLS 
z = 0.47 SN la rate of 0.42 xlO-^hfg SNe yr"* Mpc'^, 
and assume that all of this rate is generated from the 
prompt component, with a star-fo rmation density at 
z = .47 of 0.043 Mq yr"! Mpc'^ (iHopkins fc BeacomI 
l200(Tj) . this implies an upper limit for B of B < 1.0 x 10"^ 
SNe yr-i (Mq yr~i)~* , at the lower end of the SB05 as- 
sumed values. 

We can also estimate the predicted rate of SNe la in 
clusters from ellipticals and compare to the observed 
cluster rate as a further consistency check. As instan- 
taneous star-formation in massi ve clusters can b e sup- 
pressed relative to field galaxies ( Cou ch et alJ20 01^ , clus- 
ter SNe la are likely to be dominated by events from 
the old component. iGal-Yam. Maoz. fc SharonI l)2002f) 
measure a cluster SN la rate of 0.392h?o SNu". Assum- 
ing a typical elliptical mass-to-light ratio in B-band of 
5-10, this converts to a SN la rate in clusters of 3.9- 
7.8 SNe yr^^ Mq* , in excellent agreement with our A 
parameter determination. 

We also check that the A and B values derived here 
are consistent with me asured volumetric SN la rat es. We 
use the cosmic SFH of Hop kins fc BeacomI l|200fif) which 
gives the SFR of the universe as a function of redshift. 
Integrating this SFR gives the mass as a function of red- 
shift, which in conjunction with our A and B values will 
predict the SN la rate. Care must be taken that a similar 
definition of mass is used here as is used when we derive 
A i.e. the mass of stars that have died must be sub- 
tracted from the integral (see Figure [3 for the difference 
this correction makes in the mass ev olution of the cosmic 
SFH). To do this, we use the SFH of lHopkins fc BeacomI 
((2OO6) as an input to PEGASE.2 to calculate the total 
mass in stars as a function of redshift. We then use the 

14 "SNu" is one SN per century per 10^" stellar _B-band lumi- 
nosities 



A and B values from above to calculate the volumetric 
SN la rate fFigurelTHI. The relative contributions of the 
two components evolves strongly with time. The young 
component provides ~10-20% of all SNe la. at z — 0, ris- 
ing to 85% at z — 2; the exact ratios depend strongly 
on the assumed cosmi c SFH. Qualitatiyely sim ilar trends 
are also predicted bv iMannucci et al] l)2006j) . based on 
completely different analysis and fitting techniques. This 
first fraction is quite low; we note that assuming that the 
likelihood of the old component decreases with age rather 
than remains constant as in eqn. would result in a 
larger fraction of SNe being generated from the prompt 
component at 2; = 0, perhaps illustrating a limitation of 
the "A-HB" model. 

The agreement with the volumetric rates is remark- 
ably good, gi ven the S FH used in the calculation 
(^opkins & Bc acomll2006l) is completely independent of 
the derivation of A and B. All of the published rates are 
statistically consistent at the 2-tT level with the simple 
"prompt-fold" model parameters derived here. The most 
deviant point is the ~ 1 rate from Dahlen et al. (2004), 
but even this value differs from the model at only 1.8 cr. 
The A and B model used here predicts a shallow redshift 
evolution of the SN la rate, one that does not evolve as 
fast as the cosmic SFH. Of course, the distribution of 
the two types will vary considerably with local galaxy 
density. In low-mass star- forming dwarf field galaxies, 
for example, the prompt component will be the source 
for essentially all SNe la. Yet, a substantial fraction of 
SNe la will always occur in E/SO galaxies (and perhaps 
the bulges of spiral galaxies) and in clusters due the large 
amount of old stellar mass locked up in these systems. 

6. SUPERNOVA PROPERTIES AS A FUNCTION OF 
STELLAR POPULATIONS 

The results of § (S) show that the rate of SNe la in 
a given galaxy depends on both the evolved (old) stellar 
mass of that system, and on the young stellar content, pa- 
rameterized via the SFR. Indeed, a simple mass-scaling 

- assuming that the SN la rate is a simple function of 
galaxy stellar mass - is unable to reproduce the trends 
that we see. These results could have important im- 
plications for the use of SNe la as precise cosmological 
probes. At the very least, as discussed by SB05 and 
IMannucci et alJ l)200(Jl . there must be a broad range of 
delay times between a binary system formation and SN 
explosion. Whether this is interpreted as being due to a 
single channel to a SN la (i.e. perhaps single degenerate 
or double degenerate), or via a two channel path to a 
SN la (i.e. perhaps single degenerate and double degen- 
erate), a distinct possibility is that the average properties 
of these "prompt" and "old" SN la populations could dif- 
fer. In this section, we search for SN la properties that 
correlate with environment, or more particularly star- 
formation activity (or the mean age of the stellar popu- 
lation) in their host galaxies. 

One key observable affecting the utility of SNe la as 
cosmological probes is the light-curve shape/luminosity 
relationship l)Phillips|ll993D . which corrects the luminos- 
ity of SNe la according to the width of their light-curves 

- the so called "brighter-slower" relationship. In this 
paper, we parameterize t he SN la light-curves u sing the 
stretch parameter s (e.g. iPerlmutter et al.lll997|) . which 
linearly stretches or contracts the time axis of a tern- 
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plate SN la light curve around the time of maximum 
light to best-fit the observed light curve of the SN being 
fit. At low redshift, a trend of SN la absolute luminos- 
ity - or equivalently, light-curve width cast in terms of 
ATOi5(i?), the decline in i?-band magnitudes 15 days af- 
ter maximum light - versus galaxy morpho logical type 
has b een observed jiHamuv et al.lll995lll996l:lR icss ct al. 
IT999t iHamuv et al.ll200d. see lGallagher et al.H2005, for a 
recent compilation). These authors show that brighter 
SNe la (or more precisely high-stretch/low-ATOi5(i?)) 
are preferentially located in late-type (younger) galax- 
ies. 

Here, we compare the low-redshift distribution to 
that seen at high redshift in the current SNLS sam- 
ple (Figure lll() . We use the specific SFR classification 
of the SNLS host galaxies at high redshift, and take 
classifications for the low-redshift SNe hosts from ei - 
ther the literature (i Riess et al.lll999HHamuv et alJl2000t 
Ivan den Bergh et a l. 2002) or from the NASA/IPAC Ex- 
tragalactic Database (NED), simplifying onto a ellipti- 
cal (E/SO), early- type (Sa to Sbc) and late- type (Sc and 
later) classification scheme. We refit the published light- 
curves for the low-redshift SNe la with the same method 
as the high-redshift SNe la to ensure a consistent defini- 
tion of "stretch" between low and high redshift. 

The trend in Figure 1111 across galaxy type is clear 
- passive galaxies with a small (or zero) specific SFR 
tend to host lower stretch SNe. Though the low-z trend 
was already well-known, this is t he first time the ef - 
fect has been seen at high-z (e.g. iSullivan et al.ll2?)?)^ . 
presumably because the SNLS probes a wider range in 
stretch than previous high-z surveys, particular at lower 
stretches (fainter SNe). Kolmogorov-Smirnov (K-S) tests 
can be used to give some statistical guide as to whether 
these various stretch distributions are similar; these are 
summarized in Table 0] We show that at high redshift, 
the stretch distributions between SNe in passive galaxies 
versus those in both types of star-forming galaxies are 
different with >98% confidence. We also show that the 
stretch distributions for spirals at low-z and star-forming 
galaxies at high-z are essentially identical. 

The distributions in ellipticals at low-z versus passive 
galaxies at high-z are formally different at the 90% level, 
though this is caused entirely by the lack of very low 
stretch high-z SNe (s ^ 0.7). One possibility is that the 
delay times for the lowest stretch SNe la a re long enou gh 
to preclude their existence at high redshift (Howell 200 1|). 
Another is that as these SNe are the fainter subsample, 
this absence could be attributed to selection effects. 

In general, selection effects are unlikely to be respon- 
sible for the other trends that we observe. The main 
trend concerns an absence of high-stretch (i.e. brighter) 
SNe la in passive systems at high redshift; a selection or 
Malmquist bias in passive systems would lead to the op- 
posite trend, i.e. a decrease in the lower-stretch (fainter) 
SNe la in these systems. Furthermore, the trends re- 
main unchanged if the sample is culled at z = 0.6 (where 
Malmquist effects will be smaller) rather than z — 0.75. 

These trends seem highly suggestive of differing prop- 
erties between prompt and older SNe la. The passive 
galaxies, with a SFR=0 and an average older stellar 
population, host events with average stretches differing 
from those prompt SNe la found in galaxies comprised of 
younger stellar populations. Of course, galaxies with an 



average young stellar population are still likely to host 
low-stretch SNe la due to the older stellar populations 
which must be present in all but the youngest starburst 
galaxies, though the reverse is less likely to be true. We 
can attempt to identify this low stretch population in 
Figure 1111 using a simple mass-scaling technique in the 
context of the single degenerate progenitor system sce- 
nario. 

We sum the total stellar mass in all the field passive 
galaxies and in all the star-forming galaxies. We then 
scale the stretch distribution in the old galaxies by the ra- 
tio of these two total masses, and overplot on the stretch 
distributions of SNe la in the younger hosts. Subtract- 
ing this scaled population leaves just the residual stretch 
distribution remaining (Figure Clearly, other ratios 
could be used to scale the passive distribution prior to 
subtraction. We found that the residual stretch distribu- 
tion was similar when using various different mass ranges 
(e.g., 0.8 Mq to 3.0 M0) to calculate the scaling ratios. 

The resulting subtraction is remarkably clean given the 
simple assumption that was made when scaling the low- 
stretch SN distribution. The low-stretch distribution in 
passive galaxies seems able to reproduce the distribution 
of low-stretch SNe la in star-forming galaxies, within the 
statistical errors. This seems to provide supporting ev- 
idence that not only is age a key parameter driving the 
SN la rate (see §13, but that it also provides a physical 
parameter partially controlling the observed SN to SN 
stretch variation. 

7. CONCLUSIONS 

In this paper, we have examined the rates and prop- 
erties of high-redshift Type la supernovae (SNe la) dis- 
covered via the Supernova Legacy Survey (SNLS) as a 
function of the stellar mass and star-formation in their 
host galaxies. Our principal findings are: 

1. The SN la rate per unit stellar mass is a strong 
function of host galaxy specific star-formation rate. 
More strongly star-forming (later-type) galaxies 
host around 10 times as many SNe la per unit mass 
than do passive galaxies with a zero star-formation 
rate, si milar to trends observe d in the local Uni- 
verse by iMannucci et alJ l)2005D . 

2. The number of SNe la per galaxy in passive galax- 
ies closely tracks the stellar mass of the host system 
with a linear relationship. Though a relationship 
between the number of SNe per galaxy and galaxy 
mass is also seen in star- forming galaxies, it is not a 
simple linear dependence (n =~ 0.7±0.08), with an 
excess of SNe la in low-mass star-forming galaxies 
when compared to that in non-star-forming galax- 
ies. 

3. We find a clear relationship between the number 
of SNe la per galaxy and the galaxy mean star- 
formation rate averaged over the last 0.5 Gyr, even 
after removing any signal derived from the SN/host 
mass relationship found above. More strongly star- 
forming systems host more SNe la per galaxy than 
lower star-formation rate systems, with a best- 
fitting slope of rtsFR — 0.98^0;^^ 

4. By approximating the SN la rate as a bivariate 
linear function of host galaxy stellar mass (M) 
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and host galaxy mean SFR averaged over the 
last 0.5 Gyr (M„eiu), we find that the SN la rate 
in a galaxy is well represented by SNRia(i) = 
AM{t) + BMriewit), with A = 5.3± 1.1 x IQ-^"^ 
{Ho/7Qf SNe yr"! Mq^ , and S = 3.9 ± 0.7 x 10"^ 
(ifo/70)2SNe yr-i(Moyr-i)-i . 

5. We demonstrate, for the first time at high redshift, 
a relationship between star-formation in a host 
galaxy and the SN light-curve width (stretch) . We 
find that high-stretch (brighter) SNe la are exclu- 
sively hosted by star-forming galaxies, while non- 
star-forming galaxies only host low-stretch SNe la. 
We show that the SN stretch distributions in low 
redshift and high redshift spirals are statistically 
identical at 95% confidence, with the distribution 
in ellipticals identical at low and high redshift ex- 
cept for the lowest-stretch SNe la which are not 
present at high redshift. 

6. We show that the total SN la stretch distribution in 
high-redshift star-forming galaxies can be well rep- 
resented by a combination of a high-stretch compo- 
nent and a low-stretch component equivalent to the 
low-stretch distribution in passive galaxies, scaled 
by the ratio of the total mass in passive galaxies 
to that in star-forming galaxies. The indication is 
that not only can SNe la be generated from both 
old and young progenitor systems, but there is a 
systematic difference in the mean light-curve prop- 
erties of the two components. 

These conclusions could have implications for the use 
of SNe la to determine cosmological parameters. As the 
cosmic star-formation rate density shows sharp evolu- 
tion as a function of redshift, the relative mix of the two 
SN la components will change correspondingly (see Fig- 
ure IIUI) , assuming the efficiency of generating a SN la 
from a given progenitor scenario is invariant with red- 
shift. The prompt component should supply the domi- 
nant fraction of observed SNe la at high redshift, with 
the old component producing a larger fraction at low red- 
shift. The cross-over point, where the contributions from 
the two comp onents is the same, i s arou nd z ~ 0.5 — 0.9 
assuming the iHopkins fc Beacoml 1)20061) star-formation 
history. 

There are at least two direct implications for the use of 
SNe la as cosmological probes. The first is that SNe la 
should be found even at very high redshift (z ~ 3 — 4 
and above). While SN la models with significant delay- 
times would preclude the existence of SN la when the 
Universe was so young (< ~ 1 — 2 Gyr), the delay-time 



for the prompt component of SNe la is very short. The 
implication is that SNe la could therefore be used as cos- 
mological probes up until the highest redshifts at which 
stars are being formed. SNe la at z = 4 — 5 would be 
within easy range of future facilities such as the James 
Webb Space Telescope (JWST) and the various proposed 
30-lOOm class ground-based telescopes. 

The second implication is that an excellent under- 
standing of the light-curve shape/luminosity correction 
in different environments will be essential to fully exploit 
SNe la in measurements of < w >, and, in particular, 
for the more sensitive task of measuring any variation of 
w with redshift. The study and classification of SN la 
environment, using, for example, similar techniques to 
those presented here, may become as important a part 
of determining future cosmological constraints as mea- 
suring the light-curves of the SNe la themselves. Sur- 
veys which routinely obtain detailed information on the 
environment of all confirmed SNe detected will be ideally 
placed to perform studies of this nature. 
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TABLE 1 

Input and recovered masses and SFRs for our simulated galaxies 



Redshift known Redshift free 

Quantity Offset^ St. dev. Fraction of outliers Offset St. dev. Fraction of outliers 



Mass 


-0.11 


0.31 


0.042 


-0.02 


0.29 


0.101 


SFR(0.05Gyr) 


0.11 


0.51 


0.003 


0.39 


0.51 


0.032 


SFR(O.lOGyr) 


0.04 


0.47 


0.003 


0.32 


0.55 


0.026 


SFR(0.20Gyr) 


-0.04 


0.43 


0.006 


0.23 


0.55 


0.023 


SFR(0.50Gyr) 


-0.10 


0.46 


0.003 


0.11 


0.59 


0.026 


SFR(l.OOGyr) 


-0.14 


0.52 


0.031 


0.02 


0.58 


0.076 



^AU numerical values are in dex units 



TABLE 2 

Breakdown of candidate numbers and efficiencies 



Field" 


Area 


Nspec 


c 

Nphot 




Efficiency 






sq. deg. 






0.2-0.4 


0.4-0.6 


0.6-0.75 


Dl 


1.024 


25 


5 


0.314 


0.297 


0.241 


D2 


1.026 


27 


4 


0.223 


0.218 


0.181 


D3 


1.029 


28 


11 


0.331 


0.310 


0.248 


D4 


1.027 


20 


4 


0.307 


0.317 


0.274 



"Each field comprises SNe from two years of observing 
''Spectroscopically confirmed SNe la with a confidence index from 

3 to 5 iHowell et al. 2005) 
'^SNe la with either a spectroscopic r edsh ift but no definitive type, 

or solely a photometric redshift (see 5 14.21 
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TABLE 3 

Breakdown of photometric typing culls 
Criterium Nsncc 



Number over 0.2<z<0.75 116 286 

Number rejected on light-curve coverage 16 80 

Number rejected on Zspec ' ' ' 2 

Number rejected on stretch ■ • ■ 70 

Number rejected on x'^ ■ ■ ■ 75 

Number rejected on g' dispersion ■ ■ ■ 35 

Number remaining 100 24 



TABLE 4 

K-S TESTS ON THE STRETCH DISTRIBUTIONS AT LOW AND HIGH REDSHIFT 



K-S test sample 


Probability that distribution 




is the same 


Low-2 E/SO — > High-z passive 


10% 


Low-2 spirals High-z SFR>0 


95% 


High- 2: passive High-z weak SF-ing 


2% 


High-2 passive High-z strong SF-ing 


0.2% 


High-z passive High-z all SF-ing 


0.9% 


High-z weak SF-ing — > High-z strong SF-ing 


53% 



''The weak and strong star-forming galaxies are divided based on their specific 
star- formation rate; see §|F] 
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SNLS-04D1sk 
z=0.66 



SNLS-05D2hc 
z=0.36 



4" 

I 1 





SNLS-03D4dy ^ 
z=0.60 (j 



SNLS-03D4cj 
z=0.27 



+ 















Fig. 1. — Examples of the SN la host galaxy identification technique used in this paper (see S l^.ijl . Four SNe are shown. In eac h pan el, 
the SN position is marked with a cross, and each candidate host, as detected by SExtractor, has the bR ellipse over-plotted (see ? 12. 31 for 
the definition of R). The nearest host in terms of this R parameter is considered to be the correct host; SNe with no hosts inside bR are 
considered "hostless". Top Left: SNLS-04Dlsk, a straightforward case where the identification is unambiguous. Top Right: SNLS-05D2hc, 
a case where the nearest host in terms of arcseconds is probably not the correct identification. Bottom Left: SNLS-03D4dy, a case where 
no potential host is found within several arcseconds of the SN position. Bottom Right: SNLS-05D2hc, where all candidate hosts lie at 
R > 5. 
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SN Spectroscopic redshift SN Spectroscopic redsiiift 

Fig. 2. — The comparison of spectroscopic and host galaxy photometric redshift estimates for the SNLS SN la host galaxies (left) and 
the comparison between spectroscopic and SN photometric redshift estimates for all SNLS confirmed SNe (right). The solid line shows a 
1:1 agreement between Zspcc and ^photi the dashed lines show a 10% uncertainty in (l+Zspcc)- The right-hand plot, the dotted line shows 
the region above which SN e were ex c luded as SNe la during the SN completeness correction stage (S I4.2I . The SN la photometric redshifts 
follow the method of iSullivan et al.l 120061) . For the host galaxy photometric redshift, the median of Zspcc — Zphot is 0.02 with a standard 
deviation of 0.15. For the SN photometric redshifts, these values are -0.007 and 0.088. Note that the SN la photometric redshift estimates 
include a cosmological prior that precludes their use for determining the cosmological parameters. 
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Fig. 3. — The comparison of mass (M) and star-formation rate (SFR) estimates for the SNLS SN la host galaxies when the redshift 
is fixed (spectroscopic redshift; y-axes), and when the redshift is determined using the photometric redshift code, Z-PEG (photometric 
redshift; x-axes). Zero SFR systems are shown with log(SFR) = —2.5 plus a random offset for clarity. The average properties of the sample 
as determined when the redshift is known versus the photometric redshift estimate are very similar; the median difference AM is 0.026 dex 
(90% of AM lie within 0.2 dex), and the difference ASFR is 0.034 dex (90% of ASFR lie within 0.38 dex). 
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0.6 0.8 1.0 1.2 1.4 5 10 15 20 -10 -5 5 10 

Photometric stretch of LC fit g' dispersion 



Fig. 4. — The range of fitted stretch (left), hght-curve per de gree of freedom (center), and early-time g' dispersion (right) obtained 
when running the photometric redshift code of lSuUivan et all 120061) on all spectroscopically confirmed SNLS SNe la (filled histogram) and 
SNLS core-collapse SNe (open histogram) over 0.2 < z < 0.75. Cuts of s > 1.35, > 10, and g' mean dispersion >2 (see equation O*) 
are effective at removing core-collapse SNe while retaining SNe la; these cuts are shown as vertical dotted lines in the figure. 
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6 7 8 9 10 11 12 

LOG stellar mass (M®) 

Fig. 5. — The distribution of the SN la host galaxies in the SFR-mass plane. Each galaxy is coded according to its assigned type. 
Passive galaxies are shown as circles (red), normal star- forming galaxies as squares (green), and vigorous star-formers as triangles (blue). 
The black diagonal dotted line shows the division in specific star-formation rate used to sub-divide those hosts that are star-forming. The 
passive galaxies (which have a zero SFR in our models) are assigned a random SFR centered on 0.005Mq yr^^ for illustration purposes. 
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-12 -11 -10 -9 -8 

LOG Specific SFR (Moyr"* per unit stellar mass) 

Fig. 6. — The number of SNe la per unit stellar mass as a function of the star-formation rate (SFR) per unit stellar mass of the host 
galaxy. Blue points represent SNLS data-points in star-forming galaxies. The red hashed area shows the number per unit stellar mass 
as measured in the SNLS passive galaxies (assigned zero SFR in our models). Shown for comparison is the evolution in SN la rate to 
later-type galaxies observed locally by Mannucci et al. (2005), normalized to the SNLS rate in passive galaxies. The horizontal positioning 
of the Mannucci et al. data points are subject to a uncertainty when converting their galaxy types into specific SFRs. The vertical dotted 
lines show the division we use to classify the host galaxies into different types. 
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Fig. 7. — The mass evolution, as calculated by PEGASE.2, for 3 different galaxy scenarios plus the cosmic SFH of lHopkins fc BeacomI 
120061) . The different lines show the total mass in stars excluding compact objects (dotted line), the total mass in stars plus the mass in 
white dwarfs (dashed line), and the mass as calculated by simply integrating the SFH of that scenario (solid line). 
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LOG (stellar mass) 

Fig. 8. — The number of SNe la per host galaxy as a function of host galaxy stellar mass. The three colors denote SNe la in the three 
different types of host galaxy as partitioned by their specific star-formation rate. The best-fitting lines and slopes to each distribution are 
shown. For the passive hosts, a line of slope unity is also shown (dotted line). 
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Fig. 9. — The number of SNe la per host galaxy as a function of host galaxy mean star-formation rate measured over a O.SGyr period 
(see text for details). The open blue circles show the raw rate in star-forming galaxies, with the red filled squares showing the expected rate 
derived from the total stellar mass of field galaxies in each bin. The filled blue circles shows the number per galaxy after the component 
from the stellar mass is removed. The solid line shows the best-fit; the dotted line has a slope of unity. 
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1.5 
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Fig. 10. — The predicted volumetric SN la rate as a function of redshift based on the A and B values from the bivariate fit of S 15.31 
The dotted lines denotes the A "old" component (plus limits from the statistical errors), and the dashed line the B "p rompt" component 
The solid line shows the sum of the two. The filled circle is the z = 0.47 SNLS determination of the SN la rate of 'Neill et al. (2006 
Open c ircles represen t other SN la ra te determinat ions from CaoDellaro et al. (1999), Hardin et al. (2000), Pain et al. (2002), Tonrv ct a] 
120031) . TUahlen et alj J2004) and .Blanc et al.l 120041) . We conservativ ely show the statisti c al and systematic error- bars added in quadrature 
where both are given in these papers. The star-formation history of lHopkins Ik BeacomI 12009) is assumed. 
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stretch Stretch 

Fig. 11. — The distribution of the SN la light-curve shape parameter "stretch", separated according to the specific star-formation rate of 
the host galaxy. The left-hand plot shows low-redshift SNe la, the right-hand plot shows high-redshift SNLS SNe la. The typical precision on 
the stretch measure is ±0.01 — 0.02, i.e. smaller than the bin width of the histograms. LEFT: The top panel shows galaxies morphologically 
classified as spirals, the lower panel shows those SNe in elliptical or SO galaxies. RIGHT: The top panel shows galaxies with a specific SFR 
(sSFR) of log(sSFR) > -9.5, the middle panel galaxies with -12.0 < log(sSFR) < -9.5, and the lower panel log(sSFR) < -12.0. The 
vertical lines show the positions of the median stretch in each histogram. 
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Stretch 

Fig. 12. — The effect of subtracting the SN la stretch distribution in passive galaxies from the distribution in star-forming galaxies. 
The left-hand panel shows the distribution in all star-forming galaxies, the middle panel shows the passive galaxy distribution over-plotted 
scaled by the ratio of the total mass in passive galaxies to the total mass in star-formers. The right-hand panel shows the remaining 
distribution after subtracting the scaled passive galaxy distribution. 



